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A class of quasi two and three dimensional quantum lattice spin models with nearest 
(--. and next nearest neighbour interactions is proposed. The basic idea of construction is to 

introduce interactions in an array of XXZ spin chains through twisting transformation. 
The models belong to quantum integrable systems allowing explicit i?-matrix solution. 
The eigenvalue problem can be solved exactly using some symmetry of the models. 

\ 

Quantum integrable coupled spin chains with XXZ as well as Hubbard interactions were 
introduced through twisting transformation in IJ . Following that exact eigenvalue solutions of 
such models and their extensions were obtained in and ||. We generalise here this concept 
by coupling neighbouring pairs of anisotropic XXZ spin-| chains ||, which results to a family 
of quasi two (2d) and three dimensional (3d) quantum spin models with nearest neighbour 
(NN) and next NN interactions. Since the twisting transformation preserves integrability, all 
the 2d and 3d systems thus constructed turn out to be quantum integrable. Generalising 
further the useful technique of || for representing such interactions as the operator dependent 
unitary transformations, we can solve exactly the eigenvalue problem of the higher dimen- 
sional spin models. One should mention here that an integrable 2d quantum spin model, its 
corresponding classical statistical system as well as models with internal degrees of freedom 
were also obtained earlier following a different rout ||. However the twisting transformation 
due to its nice symmetry makes our approach much simpler, allowing easy construction of 
quasi two, three and in principle any arbitrary dimensional integrable spin models. 
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Let us start from spin-| operators ofi m ), a = 1, 2, 3; % = 1, 2, . . . , N; m = 1, 2, . . . , M; 
given at site (i, m) in a 2d N x M lattice with the commutation relations 

The simplest 2d quantum model we propose may be given by the Hamiltonian 

N M 

H = X! X! a {j,m) a (j+l,m) {PmO + KPml a fj,m+l) ~ Pml a fj+l,m-l)) + Pm2 (^fj im +l)^lj+l,m-i))) 
j=l m=l 

+ A m o^ )m )0-^ +1)m ) + /i.e., (2) 

where the parameters p ma involve only two independent coupling constants Q m -\m and 8 mm+ i 
between the neighbouring chains as 

Pmo = cos 2^ m _i m cos 29 mm+ i , p ml = cos 29 mT i m sin 26 mm ±\ , p m 2 = sin 26 l m _i m sin 2# mm+ i . 

(3) 

We see from the above Hamiltonian that it is asymmetric in x and y directions, which gives the 
model a quasi 2d structure. There are XX Z type NN spin interactions along the x-direction 
indicated by the sites j, while varieties of different interactions along the y direction denoted 
by the index m. 

The terms with coefficients p^i represent NN interaction involving three spins, while that 
with p m2 stands for the next to NN coupling involving four spin operators. Note that the 
anisotropy parameters are different for different chains and moreover depend on the neigh- 

pmO 

bouring chain couplings mm+ i,9 m -i m . Therefore A m — ► 1 limit can not recover the isotropic 
case. 

We intend to prove that the model represents an exactly integrable quantum system ex- 
hibiting a hierarchy of commuting conserved integrals with the Hamiltonian being one of them. 
For this we first find the associated quantum f?-matrix satisfying the Yang-Baxter equation 
(YBE) 

R<a><b>(^ — P)R<a><j>W R<bXi>(p) = R<bXj>(p) R<aXj>(ty R<aXb> 

(4) 

considering i2° 0><J - > (A) = Ylm=i R(am)(j m)(A> Vm) as the f?-matrix of M number of integrable 
noninteracting XX Z spin chains. Rfa^fj^i^Vm) corresponds to the m-th chain acting 
nontrivially on the space (h ® • • • ® I m _i ® V a <g> • • • ® hi) <8> (h ® • • • ® I m -i ® Vj ® ' • • ® hi) 
and given by its well known form || 

R(a" m )(j,rn)(^Vm) = W ( A, TJ m )I m ®I m +W 3 ( A, Vm)tfa, m ) <^j, m) + w (Vm) {° J,m) %'^) + ^m) ) ' 

(5) 
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where 

^0 + ^3 = sin(A + rj m ), w — w 3 = sin A, w = sin^ m . (6) 
We perform twisting || plus similarity transformations to construct a new i?-matrix 

R<a><j>W = F<a><j>G<a><j>R< a ><j>WG < a >< j > F <a>< j > (7) 

where 

F<aXj> = il-^(a,m)(j,m+l)j G<a><j> = JJ C(o,m)(j,m+l) (&) 
m m 

with 

F( a ,m)(j,m+1) = f(a,m)(j,m+l)(f(a,m+l)(j,m)) anC ^ C(a,m)(j,m+1) = f(a,m)(a,m+l)f(j,m)(j,m+l) (9) 

and the explicit expression through the spin operators as 

f(a,m)(j,m+l) = ex P [^mm+l°fa,m) a fi,m+l)] ■ (10) 

and with 6> mm +i replaced by a mm+ i for m )(j m+ i)- To show the new i?-matrix (|7j) to be 
a solution of the YBE (|4]) we insert it in the equation and notice that due to the specially 
designed forms @ and (|10|) of F and G, the twisting matrices can be pushed through all the 
R 0, s without spoiling their structures and canceled from both the sides. As a result we are 
left with the R° matrices only, which satisfy the YBE. For example, for shifting the F factors 
through R° in the term R <a>< b>F< a ><j>F < b><j> appearing in the equation, we notice that 
the only term in R° which does not commute with the F factors is ITm=i a ta,m) a Cb,m)- However 
using the obvious relation 

<7(a,m) exp[z0X(j, m + l)a ( 3 ajm) ] = exp[±i9X(j, m + 1)] exp[i6X(J, m + l)crf a>m) }af a>m) (11) 

with arbitrary operator X(j, m + 1), we find that the extra factor 

exp[i(^ m _i m f7y m _ 1 ) — ^mm+l (T y,OT+l))] 

produced due to the transition of i^mlfcm+il^ajm-iXj.m) is canceled exactly by the factor 
created due to F^mj^m+ij^m-ij^m). Thus R% a>< b > remains unchanged after taking the F 
factors through it. Similar reasoning holds for i?° a><J> and i2° 6>< >. The factors related to 
the similarity transformation: G7^ m ^ b m+ i)G ! (a,m)0' ) m+i) e ^c on the other hand, are partially 
canceled among themselves and the remaining ones commute trivially with the R° matrices. 
This shows that the transformed i?-matrix (0) is also a solution of the YBE, which proves 
the quantum integrability of the system and guarantees that the transfer matrix r(A) = 
(nf=i r< o><j>(A)) would generate mutually commuting set of conserved operators 
C n = (^logr(A))| A=0 1,1- 
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We construct now the Hamiltonian in the explicit form supposing H = G\ — t'(0)t~ 1 (0) 
and using the definition of r(A) along with the expressions (0) and (||). Notice an important 
property of the i?-matrix (0): 

R<a><j>(0) — cF <a>< j > G <a>< j > P <a>< j > G <a>< j > F <a>< j > 

— °F<a><j>G <a><j>G < j ><a >F<j><a>P<aXj> = c P<aXj>i (12) 

rM 



which follows easily from that of i?° a>< > (0) = 11^=1 s ™ VmP(a,m)(j,m) — °P<axj> an d the 
symmetry F <j><a> = F< a 1 ><i> , G <j><a> = G <a><j> . Using this along with the prop- 
erty of the permutation operator like P<a><j>R'< a xj+i>(®) = P'<jxj+i>(fyP<a><j>, an d 
P <axj> = 1 we § et 

M 



and 



r (°) = c ( ( II P (j,m)U+i,m))P<jxj+2> ■■■) tr a (P <a><j> ) (13) 

V m=l J 



N M m 

r '(°) = C XXX II ( F UMj+hl+l) G Um+hl+l)) R tiZ)ti+l,rn)(Q) II P(j,n)(j+l,n) 
j=lm=l l=m—l n^m 
m 

II P(j,iKj+^Hi) G (j,i)(j+i,i+i))P<j><j+'2> • • • x tr a (P <a xj>), (14) 

l=m— 1 

yielding 

iV M m 

TT \ \ am jdIxxz (r\\ JD (om \ — 1 ™ TT 771 ri 

n - 2^ 2^ ^jj+l K (j,rn)(j+l,m){v)^j,m)(j+l, m ){Jjj + l) , ~ H ^ <J,l)(j+l,l+l) ^(j.OO'+l^+l) • 

j=l m=l Z=m— 1 

(15) 

Expressions (|5|) , fllOD through spin operators and their properties ([□]) reduce (|T5| ) finally to a 



family of 2d quantum spin models 

N M 

^ = EE a {j,rn) a (j+l,m) eX P(*[(^mm+l + a ram+l)°fj, TO +l) _ (#m-lm + am-lm)^ + l, m -l) 
j'=l m=l 

+ {Omm+1 — a mra+l)°Yj+l,rn+l) ~~ ($m-lm ~ a m~lm)C( 3 j >m _l)]) + ^•rnC( ; / jTO )0'(j+l im ) + ^- c -i 

(16) 

parametrised by the set {0 mm+ i}, {a mm+ i}. This also establishes that the family of models fllTf ) 
belongs to the hierarchy of quantum integrable systems. Choosing the coupling parameters 
9, a differently we can generate quasi 2d integrable models with rich varieties of interchain 
interactions. The simplest choice a mm+ i = 9 mm+ i, using (c"p am )) 2 = I clearly yields the 
Hamiltonian (H), while au + \ = (—l) l 9u +1 leads to the models like 

N M 

ff = EE ( a u, m )%+hm) S m + A m af lm) af J+1<m) ) + /i.e., (17) 

j=l m=l 
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with S m taking alternate expressions as 

S m = ex P(2«[6'mm+lO'fj im+ l) ~ ^m-lm^m-l)])) 

S m +i = exp(2i[6 m+ i m+2 crfj +ltm+ 2) — &mm+l^{j+l, m )])- (18) 

One may also choose further the M number of parameters 6 mm+ i differently to get other types 
of interacting 2d models. 

Next we address the eigenvalue problem: H | $ >= E | $ > and show that by exploiting 
the symmetry of the operators E (j !m )(j+i,m+i), (% m )y+i, m +i) appearing in (0) and using the 
known results for the XXZ spin chains one can solve this problem exactly. Indeed, using the 
relations like (|TTJ) one can conclude interestingly that our Hamiltonian ( |T5| ) can be reduced to 
the form H = SHqS^ 1 , where H is a collection of exactly solvable XX Z chains and S is an 
unitary operator given by 

,m)(fc,m+l) J (i,m)(j,m+l)\J (i,m+l)(j',m) (19) 

k,m,i,j(i<j) 

with f d ,f a as in fllPp. Note that for the particular m = 2 case and a = the result is 
consistent with B]. Therefore the eigenvalue problem of the quasi 2d Hamiltonian (0) can be 
reduced to that of H with the same energy spectrum E = Y^m YJk™ C0S Pk m , where M m are 
the number of spin excitations in the m-th chain. The eigenfunction on the other hand may 
be given by S(Yl m \ >), where | > is the known solution for a single XX Z chain 

solvable by the Bethe ansatz 0. Note that, though our model can be solved through integrable 
XX Z chains, the results are not really equivalent. Indeed, the energy spectrum depends on 
the values of the unknown momentum parameters pk m , which in turn can be determined by 
the Bethe equations and these equations are not the same for the two models. In the 2d model 
due to the operator S extra factors like exp[z(6 |mm+1 M m+1 — 6 m ~ lm M m _i)} appear in the 
determining Bethe equations signifying interactions between the neighbouring chains and this 
gives the set of momentum parameters a different value. The nature of the eigenfunction is 
also changed due to the same reason. At different sectors in the configuration space different 
phase factors arise in the wave functions due to the operator nature of S and at the boundaries 
of the sectors the wave functions suffer phase jumps resulting the occurrence of discontinuity 
at coinciding points. This unusual feature was observed also in case of the twisted Hubbard 
model 0. 

We extend now this idea to construct quasi 3d quantum spin models by defining the R- 
matrix (^) with transforming operators 

p _ TT p( m ) pip) n _ TT ni m ) niv) 

r <a><j> — 1 1 r (a,m,p)(3,m+l,|)) r (a,m,p)(j,m,p+l) ) ^<a><j> ~ H ^(o.m ,p) (j,m+l,p) (a,m,p) (j ,m,p+l) 
m,p m,p 

(20) 
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Here as before the index j = 1,...,N denotes the site number in the x-direction, m = 
1, . . . , M stands for the chain number along the y-direction, while the additional index p = 
1, . . . , L indicates the layer number along the ^-direction. Therefore we may extend our above 
definitions to the 3d case as 

p(p) __ j-Opp+i l f^pp+i \— 1 

(a,m,p)(j,m,p+l) J (a,m,p)(j,m,p+l) \J (a,m,p+l)(j,m,p) ' 

and 

(-<{P) f a PP+l £ a PP+l (21) 

(a,m,p)(j,m,p+l) J (a,m,p)(a,m,p+l)J (j,m,p)(j,m,p+l) \ I 

with 

f(a%.!p)(j,m,p+l) = ex P[^PP+l a (a,m,p) a (j,m,p+l)\- (22) 

for fixed chain index m and 9 pp+ i as the interlayer coupling. Analogous relations as (|2l|) hold 
also for Ffam^u m +i P ) with fixed layer index p. All the above arguments for proving the 
i?-matrix as the YBE solution also go through for the present extension and one can generate 
the related Hamiltonian of a family of quantum integrable 3d spin models. In the simplest 
case of a a b = 9 a b the explicit form of such models may be given as 

N M L 

H = E E Ef^.m./O+l.m.p)^ 1 + A m 4,™,/(J+I,™,ri) + ftX ' ( 23 ) 
j=l m=l p=l 

with the interchain and interlayer interactions given as 

*^m,p = ex P(2^[(^mm+lC"f jr - im +i iP ) ~ ^m-lm cr f :7 +l,m-l,p)) + {@PP+l a (j,m,p+i) ~ ^p-lp a (j+l,m,p-l))\) • 

Note that due to the factorised form of 

a m,p " J (m+l,p)\' J (m-l,p)J u (m,p+l) W(m,p-1)/ 

with 

S (m, P +i) = expi^PP+itfj^p+i)]) = cos 29 pp+i + icr U,m, P +i) sin 26 pp+1 (25) 

etc. such interactions take place between NN and next NN chains in the same plane (with 
fixed p) as well as between NN and next NN layers in the same chain (with m fixed) and 
finally between such terms themselves. 

The eigenvalue problem of such 3d models can also be solved similar to the 2d case by 
transforming the Hamiltonian to the noninteracting array of XX Z spin chains through unitary 
transformation. 

Thus we have constructed and solved a class of quasi 2d and 3d quantum spin models 
and shown also their exact integrability. It should be obvious from our construction that 
the approach allows further extension of such integrable models to any arbitrary dimension. 



(24) 
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Unlike the long ranged systems 0, the present class of 2d models includes only nearest and 
next nearest neighbour interactions. Therefore such models, possibly with more interesting 
twisting operators, might be helpful in constructing integrable higher dimensional physically 
significant models, which could not be achieved through the approach of tetrahedron equation 



10] . It would also be important to construct higher dimensional Hubbard like models through 



similar technique [|TT|j . 
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